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The phase diagram of the t — J — V model is discussed using a 1/N expansion in terms of X operators. It is 
shown that a flux phase of d-wave symmetry is stabilized by the Coulomb interaction at intermediate dopings 
and competes with d-wave superconductivity. Since the flux wave instability is stronger than the superconducting 
one optimal doping is essentially determined by the onset of the flux phase at zero temperature at the doping 
S c ~ 0.14 for J/t — 0.3. For 8 < S c the flux phase coexists with superconductivity at low and exists as a single 
phase at higher temperatures. Due to the different origins of the two instabilities the flux phase boundary and 
especially S c is much less sensitive to impuritiy scattering than the boundary for superconductivity in agreement 
with experiments in Zn doped La — 214 and (Y, Ca) — 123. 



1. INTRODUCTION 

The superconducting transition temperature 
T c of many high-T c superconductors can strongly 
be suppressed by Zn impurities. Considered as 
a function of the doping 8 T c vanishes finally at 
a doping Sqcp which is the T — endpoint of 
a line separating a metallic phase at higher dop- 
ing from a pseudogap phase at lower doping in 
the T — 5 plane 0. The microscopic nature of 
the pseudogap phase is presently not clear. Since 
its gap has d-wave symmetry like the supercon- 
ducting gap it has been suggested that Cooper 
pairs without long-range phase order are respon- 
sible for this phase. On the other hand there is 
evidence for the occurrence of strong fluctuations 
in the underdoped regime due to the proximity of 
phases different from the superconducting phase. 
At zero temperatures these phases may exhibit a 
change from long- to short-range order at Sqcp 
so that Sqcp is the location of a quantum critical 
point. One aim of the present investigation is to 
study the t— J— V model at large N to see whether 
such a quantum critical point scenario applies in 
this case and what the relevant additional phase 
is besides of the superconducting one. Another 
aim is to study the relation between these two 
phases. 



2. COMPETITION BETWEEN SUPER- 
CONDUCTIVITY AND FLUX 
PHASE 

We consider a t-J-V model with N degrees 
of freedom per lattice site on a square lattice. 
Its Hamiltonian can be written in terms of X- 
operators as 



H 



Et 
1 



ij v p0 v 0p 

N Xi X > 



E 



Jjj vPI yIP 



\ " / Jij , *y \ yPP yqq 



(1) 



The internal labels p,q... consist of a spin label 
distinguishing spin up and spin down states and 
a flavor label counting N/2 identical copies of the 
original orbital. The coupling constants tij and 
Jij are confined to nearest neighbors i,j and sim- 
ply denoted by t and J, respectively. The X op- 
erators satisfy the commutation and anticommu- 
tation rules of Hubbard's X operators for all N. 
Moreover, the sum of the diagonal operators is 
equal to N/2 at each site meaning that only N/2 
out of the N states per site can be occupied simul- 
taneously. The first three terms represent the t-J 
Hamiltonian, the last term a screened Coulomb 
interaction appropriate for two dimensions and 
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taken from Ref. |J2|. In the following we express 
all energies in units of t and all lengths in units 
of the lattice constant of the square lattice. The 
strength of the Coulomb interaction will be char- 
acterized by its value between nearest neighbor 
sites V nn - 

The self-energy £ of the one-particle Green's 
function is independent of frequency at large N 
and has the general form || 

E(k,q)=E(k)iV c *(q) + 0(k,q). (2) 

S(k) denotes the self-energy in the normal state. 
<p is an additional contribution describing a new 
ground state characterized by the modulation 
wave vector q. For the determination of phase 
boundaries (\> can be considered as infinitesimally 
small. Calculating explicitly the self-energy from 
Eq.(|l|) one finds 

0(k,q)=5^/ a (q)F a (k,q), (3) 

a 

with F(k, q) = (t(k — q), 1, Jcos(fe x ), Jsin(k x ), 
Jcos(ky), Jsm(k y )), and f a satisfying the ho- 
mogenous system of equations 

Y, [<W - oa/?(q)] Mq) = o • (4) 



The matrix elements a a p(q) in Eq. (^) are de- 
fined by 

c p n 

5o(p,«w„)5o(p - q, iu n ), (5) 

with G(k, q) = (1, i(k) + F(q) - - cos(fc x ), 
- sin(fc x ), - cos(fc y ), - sin(fcj / )), and g (p, = 
l/(iu> n — E(p)). The boundary between the nor- 
mal state and a possible new phase corresponds 
to a nontrivial solution for fp of Eq.(Q). In 
this way both the modulation vector q and the 
components of the order parameter in the six- 
dimensional order parameter space is determined. 

Putting first V = and decreasing doping or 
temperature from large values we find for J < 1 
a first instability of the normal state with respect 
to the flux phase order parameter 

(f'FL (k, PQ) ~ i(cos{k x ) - cos(k y )). (6) 
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Figure 1. Incommensurate (solid line) and com- 
mensurate (dashed line) phase boundaries in the 
T — 5 plane. Insert: Instability vector Q = (l,x)ir 
as a function of T. 

The solid line in Fig. 1 shows the correspond- 
ing instability line for J = 0.3 in the T — 8 plane. 
Writing the modulation vector as Q = 7r(l, x) the 
inset of Fig. 1 shows function of temper- 

ature along the instability line. For T > 0.015 
the flux phase is commensurate with the wave 
vector Q c = tt(1, 1), for T < 0.015 one wave vec- 
tor component is incommensurate. Fixing Q to 
the commensurate value Q c for all temperatures 
the instability line is given for T < 0.015 by the 
dashed line in Fig. 1 describing a reentrant be- 
havior of the normal state between S = 0.12 and 
0.14. This feature is an artefact caused by the 
assumption of a commensurate phase at all tem- 
peratures. At zero temperature the critical dop- 
ing S c is a monotonic function of J, with S c = 
for J = 0, S c ~ 0.14 for J = 0.3 and 5 C - 0.20 for 
J = 0.6. Thus 5 C is near the experimental value 
for optimal doping and denotes a quantum criti- 
cal point because the incommensurate flux phase 
has short-range order for any finite temperature 
if fluctuations are taken into account. 

The flux phase is unstable at very low dopings 
with respect to dimers [Q. Another instability 
line, indicated by the dot-dashed curve in the up- 
per panel of Fig. 2, describes a diverging com- 
pressibility of the normal state, i.e., the region to 
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Figure 2. Phase diagram for J = 0.3 and 
different Coulomb interaction strengths V nn . 
N,FL,SC,ICDW denote the normal, flux, su- 
perconducting, and CDW phase, respectively, the 
dot-dashed line in the upper panel a diverging 
compressibility in the normal state. 

the left of this line would show phase separation 
in the normal state. Since the flux phase falls in- 
side the phase separation regime, the solid line in 
Fig. 1 has no physical meaning in the case V = 0. 

The instability line of the normal state towards 
superconductivity has been discussed in Ref ||. 
In the BCS-approximation T c is ~ exp(l/Xi), 
where Xi is the lowest eigenvalue of the static ker- 
nel of the linearized gap equation and i denotes 
one of the 5 representations of the point group 
Civ of the square lattice. Only A3 correspond- 
ing to d-wave symmetry exhibits strong negative 
values. It even diverges at 5 C due to critical fluc- 
tuations caused by the incipient flux phase. The 
numerical solution for T c , which also takes into 
account retardation, however, shows that criti- 
cal fluctuations have only little influence on T c 
and that T c is mainly determined by the instan- 
taneous term. Keeping only this term the solid 
line in the upper panel of Fig. 2 shows T c as a 
function of S. Interesting is that the supercon- 
ducting unlike the normal state is stable against 
phase separation. 

Taking also the Coulomb interaction V into ac- 



count phase separation is no longer possible but 
charge density waves (CDW) may be stabilized. 
The two lower panels of Fig. 2 show, however, 
that CDW's may be only stable at rather low dop- 
ings far away from optimal doping and are thus 
of less interest for the following. In order to be 
able to discuss fully the interplay between super- 
conductivity and the flux phase we treat the cor- 
responding order parameters in a nonlinear way, 
assuming only that the flux phase is commensu- 
rate. The self-consistent equations are then 4x4 
matrix equations which are most conveniently for- 
mulated using a Nambu representation with 4 
states. 

The middle panel of Fig. 2 shows the phase 
diagram for V nn /J = 0.5. For S > S c T c is 
monotonically increasing with decreasing 6 sim- 
ilar as in the upper panel, only its absolute value 
is reduced because of the repulsion V nn . Below 
the onset of the flux phase at S c superconductiv- 
ity and the flux phase compete with each other 
because both have d-wave symmetry. Since the 
nesting instability leading to the flux phase is 
stronger than the superconducting instability the 
flux phase suppresses strongly T c with decreasing 
8 below S c . As a result T c assumes a maximum 
just below 8 C and decreases both towards lower 
and higher dopings. Optimal doping is thus de- 
termined at large N by the onset of the flux phase 
at S c . For 5 < 5 C the flux and the supercondcut- 
ing phase coexist for T < T c whereas the flux 
phase alone is stable for T c < T < T*, forming 
there a pseudogap phase with a d-wave gap in the 
single particle spectrum. The lower panel in Fig. 
2 shows the phase digram for V nn /J = 1. The 
Coulomb repulsion is now strong enough to sup- 
press completely T c . The region of CDW states 
has been shifted further towards smaller doping 
and the flux phase covers a large region to the left 
of the normal state. 

3. PHASE DIAGRAM IN THE PRES- 
ENCE OF IMPURITY SCATTERING 

Both the flux and the superconducting phase 
have d-wave symmetry. Their phase boundaries 
T* and T c , respectively, thus should be sensitive 
to impurity scattering. Experimentally this holds 



4 



in the case of Zn impurities for T c but practically 
not for T* 0. From a theoretical point of view 
T c and T* may depend differently on impurity 
scattering because the underlying instabilities are 
quite different. For instance, nesting effects play 
an important role in the case of T* but not in the 
case of T c . 
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Figure 3. Solid lines: Suppression of T c by im- 
purity scattering, calculated with the scattering 
rates T = (squares), 2 • 1CT 3 (circles), 4 • 10~ 3 
(diamonds), and6-10~ 3 (triangles). Grey region: 
Corresponding variation of the transition temper- 
ature T* to the flux state. 

In a simple approximation, the effects of impu- 
rities in the normal state can be taken into ac- 
count by introducing a renormalized frequency 

iw n = iuj n + iF-r^T' ( 7 ) 

\UJ n \ 

where T is a scattering rate, here used as a free pa- 
rameter proportional to the impurity concentra- 
tion. Throughout the flux phase the self-energy 
due to impurity scattering is still diagonal in the 
4x4 Nambu representation because the flux order 
parameter does not couple to the impurities. The 
constant T in Eq. (]?]) accounts in a phenomeno- 
logical way both for weak and strong potential 
scattering by impurities. The interesting dop- 
ing region for superconductivity is in our model 
rather far away from the Van Hove singularity, so 
that the band can be assumed to be structureless 
and r to be independent of frequency. 



The solid lines in Fig. [3| show numerical results 
for T c as a function of doping S for different scat- 
tering rates T, using J = 0.3 and V nn = 0.5J. 
These curves illustrate the suppression of T c with 
increasing scattering rates Y = 0, 2-10 -3 , 4-10" 3 , 
and 6 • 10" 3 . The corresponding changes in T* , 
determining the phase boundary between the nor- 
mal state and the flux state, are depicted in Fig. 
H by the grey region. The chosen values for T 
correspond roughly to T ~ 1.0T C at optimal dop- 
ing, and to r ~ 1.5T C in the strongly underdoped 
region, interpolating between the weak- and the 
strong-coupling regimes. One important result 
of Fig. U is that the flux phase boundary T* is 
only slightly shifted by impurities, in spite of the 
strong suppression of the superconducting critical 
temperature. In particular, the critical doping S c 
at zero temperature is almost completely inde- 
pendent of the impurity scattering rate. Since 
in our approach the maximum of T c as a func- 
tion of doping is essentially determined by S c this 
means that the T c (5) curves shrink to S c with in- 
creasing scattering rate which is a characteristic 
feature of Fig. ||. Interpreting Fig. || in terms of 
a quantum critical point scenario means that the 
corresponding critical doping Sqcp is given by S c 
and that Sqcp is nearly completely independent 
of the impurity scattering rate. The curves in Fig. 
H are in excellent agreement with the correspond- 
ing experimental curves in Zn doped (Y,Ca)-123 
and La-214, as given in Fig. 2 of Ref. Q]. 
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